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Abstract 

The phase space of an integrable, volume-preserving map with one action and d angles is foliated 
by a one-parameter family of d-dimensional invariant tori. Perturbations of such a system may 
lead to chaotic dynamics and transport. We show that near a rank-one, resonant torus these 
mappings can be reduced to volume-preserving "standard maps." These have twist only when the 
image of the frequency map crosses the resonance curve transversely. We show that these maps 
can be approximated — using averaging theory — by the usual area-preserving twist or nontwist 
standard maps. The twist condition appropriate for the volume-preserving setting is shown to 
be distinct from the nondegeneracy condition used in (volume-preserving) KAM theory. 



1 Introduction 

Volume preserving maps are appropriate models for many systems including fluid [CFP96, MGPM99, 
ISVL011 lRML+031 ISCVH041 IATPM061 IM,TM08j and magnetic field line flows [TH851 IGre931 IBDS98] 
and even the motion of comets perturbed by a planet on an elliptical orbit [LS94J. One fundamen- 
tal question relevant for applications is to understand transport in these maps [PF881 1RKKCA93, 
ICFP941ILM09] , 

Transport is prohibited between any two regions that are separated by a codimension-one sur- 
face. For volume-preserving maps, these dividing surfaces are often tori. Even when there is no 
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dividing surface, we may imagine, as in the two-dimensional case [MMP84J, that remnants of such 
surfaces — if they exist — may impede transport. Similarly, transport may be reduced because of 
the "stickiness" of invariant tori, a phenomenon well-established numerically for area-preserving 
maps [M086, CK08, Ven09] and for nearly integrable Hamiltonian systems P93, PW 94]. if not 
yet completely understood. For three-dimensional maps, stickiness has been observed as algebraic 
decay of exit time distributions for some systems |MJM08] but not for others |SZ08j. 

Codimension-one invariant tori are a common feature of the dynamics of three-dimensional 
volume-preserving mappings. Indeed these maps could be regarded as "integrable" when their 
phase space is foliated by a family of two-dimensional tori on which the dynamics is conjugate to 



a rigid rotation. We can think of this case as a map with two angles and one action [PF88J. KAM 
theory implies that two-tori are robust features of nearly-integrable, one-action maps [CS90, Xia92]. 
By contrast, in maps with one angle and two actions (which might also be called integrable), the 
invariant circles are apparently not robust [MezOlj . 

Tori also arise naturally through bifurcation. For example, a fixed point with multipliers 
^ e 27na^ e -27uo^ generically undergoes a "saddle-center-Neimark-Sacker" bifurcation that can cre- 
ate an invariant circle surrounded by a family of two-tori |DM09] . This bifurcation has also been 
studied in the dissipative case where it can lead to attracting invariant tori as well as strange attrac- 
tors [BSV08a, BSV08bj- Tori are often destroyed by resonant bifurcations, and this phenomenon 
is the topic of the current paper. 

We begin, in ^2j with a short discussion of integrable, volume preserving maps with one action 
and d- angles, and the frequency map, $7 : R — > WL d , that characterizes their dynamics. Our first 
goal is to obtain a "standard" volume-preserving mapping that describes, as Chirikov's standard 
area-preserving mapping does, the destruction of these invariant tori as the strength of the resonant 
terms grows. We begin in ^3] with a numerical study of a simple, three-dimensional example; similar 
examples have been studied in CK()7j. As we will see, this example exhibits many of the features 
of general one-action maps, near a "rank-one" resonance. 

In ^4] we consider more generally a small perturbation of the integrable, one-action map. Zoom- 
ing in to a neighborhood of a rank-one resonance gives rise to a map that has a phase space with 
d — 1 fast angles, one slow resonant angle and one slow action. The three-dimensional version of 
this map (when there is one fast angle) is — to first order — equivalent to the example studied in 

The main tool for nearly-resonant dynamics is averaging theory. We will show, in £j5j that 
volume-preserving normal form theory can be used to transform to new variables that approximately 
decouple the slow dynamics. This results in an approximate translation symmetry along the fast 
angles. The normal form can be computed to arbitrary order in the perturbation strength and 
can be made volume-preserving. The approximate symmetry becomes an "exact" symmetry if the 
high-order, fast fluctuations are discarded, or equivalently if this system is averaged over the fast 
variables. A symmetry reduction then leads to a map on the slow phase space that is again volume 
preserving. As is usual for averaging theory, the resulting averaged system is shown to closely 
approximate the full dynamics, but only on a finite, but long, time scale, see Thm. [5] 

Nevertheless, the qualitative features of the decoupled, averaged map agree surprisingly well 
with the numerical results of £j3j provided that one is far from "resonance overlap." In particular, 
when there is one fast angle, the reduced slow map becomes Chirikov's area-preserving standard 
map, or its non-twist cousin, depending on properties of the frequency map Q. The (2D) resonance 
structure and Chirikov's (2D) overlap criterion agree well with slices of the full three-dimensional 
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dynamics. 

Violation of the twist-condition for area-preserving maps leads to meandering curves [HH84. 
Sim98, WAFM05J and twistless bifurcations [DMSOO]. We will see in ffl and more generally in 
^4] that the transversality conditions whose violation leads to the same dynamical phenomena in 
higher- dimensional, volume-preserving maps is not the violation of a Kolmogorov nondegeneracy 
condition, but instead the condition that the image of Q be transverse to the resonance surface in 
frequency space. Our condition is weaker: KAM theory may still hold at a tangency. This is an 
essential difference to the symplectic case and shows that KAM theory "knows" little about the 
qualitative dynamics. 

A similar analysis (blowing-up, decoupling, and averaging) applies to rank-two resonances, as 
we discuss in £j6j though in this case there are three slow variables. The resulting reduced, three- 
dimensional, slow map is becomes a quasi-periodic version of the standard map. 



2 One-Action Maps 

Though applications of volume-preserving dynamics mainly occur in three-dimensions, where two- 
dimensional tori are barriers, it is easy to generalize to d-tori in d + 1 dimensions. The unperturbed 
model is the integrable family of maps /o : T d x IR — ?- T d x IR. given by 

fo(x, z) = (x + Q(z) mod 1, z) . (1) 

Such a map is called a one-action map: the variable z plays the role of the action, and the d-angles, 
Xi ,i = 1, . . . d, are taken to have period one |CFP96[ ICFP991 IMezOl] . This map can be taken 
to be a model for the behavior in the neighborhood of a smooth, codimension-one invariant torus 
on which the dynamics is conjugate to rigid rotation. The frequency map (or perhaps, using the 
terminology from flows, the "frequency-ratio" map), 

n : R R d 

traces a curve in the frequency space, oj = ft(z), as the "action" varies, as sketched in Fig. [TJ If we 
fix a lift for the periodic component of the map /o then £l(z*) is the rotation number for orbits on 
the torus z = z* . 

A rotation vector oj G M rf is said to be resonant if there exists a nonzero (m, n) E Z d x Z such 
that 

m ■ oj = n ; (2) 
if there are no such integers, then oj is incommensurate. The collection of d — 1 dimensional planes 

11 = G R d : m ■ oj = n, (m,n) € Z d+1 \ {0}} 

is the resonance web; it is a dense subset of M. d . 

Our goal is to understand the effects of perturbation and resonance on the tori of Q. As we 
will see, it is important to distinguish between resonances of different ranks. For a given oj, the 
resonance module is the sublattice of "L d that corresponds to all integer vectors m that satisfy ^ , 
namely, 

C(oj) = {m£Z d :m-oj£Z}. (3) 
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Figure 1: Resonance web, thickened by the Diophantine condition Q and a frequency map f2 for a one-action map 
on I x T 2 . The complement of the web is the positive measure cantor set of Diophantine rotation vectors. 

When oo is incommensurate C{oS) is trivial: it contains only the origin. When this set is r- 
dimensional, i.e., is the integer span of r independent vectors, then the resonance has rank r. 
When Q(z*) is r-resonant, the invariant ci-torus z = z* of ([I]) consists of families of orbits that are 
dense on (collections of) (d — r)-tori, see App. |Aj 



The d-dimensional invariant tori of ([!]) are examples of rotational tori; generally, a rotational 
torus is a set that is homotopic to the horizontal torus {(x, 0) : x G T rf }{^] Upon perturbation, 
many of the rotational invariant tori of ([!]) will be destroyed. In particular, the resonant tori are 
fragile — we expect that they will be immediately destroyed and replaced by a finite number of lower 
dimensional tori and perhaps by families of nonrotational, d- dimensional tori. Conversely, KAM 
theory teaches us that for a torus to be robust it must be not only incommensurate, but sufficiently 
incommensurate — typically it must have a Diophantine frequency vector. The set of Diophantine 
frequency vectors is 



where \m\ is any norm, for example the maximum norm. Even though 1Z and T> are disjoint and 
TZ is dense in R d , V is a positive measure Cantor set for each [i > d and sufficiently small c > 0, 
see Fig. [T| 




P(c, /i) = {uj G M d : \m ■ uj — n\ > c\m 



» V(m,n) GZ d+1 \{0}} 



(4) 



1 Rotational tori need not be graphs over the angles, see §3 
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KAM theory can indeed be applied to volume-preserving perturbations of the one-action map 
([!]); it implies, under certain smoothness and nondegeneracy assumptions, that there exists a set 
of invariant, Diophantine <i-tori (with \x > d 2 + d — 1) whose measure approaches one as the 
perturbation goes to zero |CS90| IXia92j . In these KAM-type theorems a sufficient nondegeneracy 
condition is that the Wronskian of D£l must be bounded away from zero. 

The curve Q(z) will typically intersect planes corresponding to 1-resonances for a dense set of z 
values. In this paper we want to emphasize that the dynamical behavior of the perturbed map near 
these intersections depends in detail on the rank of the resonance and on whether the intersection 
is transverse or tangent. When the intersection with the resonance plane is transverse, we will 
show in S|4] that the perturbed map can be approximated by the area-preserving standard map that 
drives additional fast angles. By contrast, when the curve Q(z) is tangent to the resonance plane, 
the map can be approximated by the area-preserving standard nontwist map that drives additional 
fast angles. 

3 Example: Transverse and Tangent Resonances 

We start by studying the dynamics of an example of a perturbation of /o, a map (x' , z') = f £ (x, z) 
given by 

x = x + ft(z') mod 1 , 
z = z + eg[x) , 

with g(x + m) = g(x) for all x G T d and m G 1j d - Here we study the case d = 2, but in ^4] we 
consider arbitrary d as well as a more general perturbation of ([!]). 

The map ^ is one natural generalization of the Chirikov standard map (a one-action, one- 
angle map) and its symplectic cousin, the Froeshle map [Mci92j (a two-action, two-angle map). It 
is volume-preserving with the standard volume form Adxi A ... A dxd A dz\ it is not hard to see that 
the Jacobian, Df £ , has determinant one. The map §5§ has no rotational invariant tori unless the 
function g has zero average: 

g{x)dx\ . . . dxd = . (6) 



/ 



This is precisely the condition that f £ is an exact volume-preserving map, or equivalently that it 
has zero net flux. Since we are interested in the persistence and destruction of invariant tori, we 
will assume that the zero flux condition ^ holds. 

For the three-dimensional case, the image of the frequency map is a curve in M? that can locally 
be approximated by a parabola. As we will see more generally in £j4j coordinates can be chosen 
near a rank-one resonance so that the frequency map becomes 

fi(z) = (7 + z, -5 + /3z 2 ) . (7) 

This parabola is sketched in Fig. [2] The twist condition of the volume-preserving KAM-type 
theories in this case becomes (3 ^ 0. 
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Figure 2: Resonance lines m ■ lo = n for |m| < 3, and the image of the frequency map |7| for j3 = 2, 7 = y/5 — 2 ; 
0.236, and two values of S that give rise to tangencies of fi with the forced resonances of |9j. 



Using the frequency map ([7]), ([5]) takes the "standard" form 

x = x + z + 7 , 

y> = y-S+ (3z' 2 , 

z ' = z + eg(x,y) . 



(8) 



The forcing function g{x) in Q is periodic; therefore it can be expanded in a Fourier series. Since 
g(x) has zero mean ([6]), its (0,0) Fourier component must vanish. In this section, for simplicity 
and to illustrate the general theory of $4l we choose g to have just three terms: 



g(x,y) = — asin(27rrr) — 6sin(27ry) — csin(27r(x — y)) . 



(9) 



The terms with amplitudes a, b, c represent forcing at the resonances (1, 0, n), (0, 1, n), and (1, —1, n), 
respectively. These forced resonances, m ■ Q(z*) = n, occur at positions z* listed in Tbl. [TJ 

Since Q(z) is a parabola, it may intersect a resonance line m ■ uj = n more than once; however 
generically, it will cross transversely. When e<1 the dynamics near a transverse crossing of a rank- 
one resonance can be analyzed by expanding the map near z = z* and performing an appropriate 
average over the nonresonant, fast angle using the averaging theorem presented in ^5} 
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(m, n) Amplitude 



M 



5t 



(1,0, n) 
(0,1, n) 
(1,-1, n) 



a 
b 



n 



± 



-7 

<5+n 

/S 



J_ 

2/3 



1 ± ^/l + 4/3(5 + 7 



1 
i 

2/3z* 
1 

l-2/3z* 



none 
— n 



Table 1: Forced resonances, m ■ Q(z*) = n for ([8]), their effective masses M, and values of 5 for a 
tangency. 



For example, the (1,0,0) resonance occurs when Q x (z a ) = 7 + z a = 0. Generically the second 
frequency Q y (z a ) = f3^f 2 — 5 = loq will be irrational so that the resonance has rank one with 
resonance module 

£{(0,u )) = {{k,0):keZ} , ujq^Q- 
To expand near this resonance, define new variable^] 

(^V,C) = (y,x,e- 1/2 (z-z a )) . 

In these coordinates, (181) reduces to 



e' = e + ^o + o( v / i) 

7/ = rj + V^C , 

C' = C + v%(^) • 



(10) 



When £<1, the angle £ is rapidly rotating relative to 77, so that the forcing terms in ([9]) proportional 
to b and c rapidly fluctuate relative to the term proportional to a. Thus one would expect that 
the dynamics of the slow variables (77, Q could be approximated by simply averaging over the fast 
angle, giving the two-dimensional area-preserving map 



7/ = 7] + y^C' , 

( = ( — \/ea sin(27T7/) 



11) 



Indeed, this is verified in ^provided that the fast frequency, ojq, is Diophantine Q. The transfor- 
mations of ^5] — to first order — then give the map (11). Theorem [5] implies that an orbit of (11) is 
0(e)-conjugate to the projection of an orbit of the original map onto the slow variables (77,0 on 
the time scale t ~ C(e -1 / 2 ). This is valid so long as the orbit of ( 11 ) stays in an 0(1) neighborhood 
of C = 0. The quality of approximation can be increased to 0(e Tl / 2 ) for arbitrary n by the iterative 
procedure described in $5} 



Note that (11) is symplectic. Indeed, when e <C 1, it is approximately the time y/e map of the 



one-degree-of-freedom pendulum Hamiltonian 



%C;M,o) 



2M 



2vr 



cos(2-7r?7) 



12) 



2 We switch the order of x and y to agree with the general notation of [JH] 
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with momentum £, coordinate rj, an effective "mass" M = M a = 1 and amplitude a. 

Consequently, in a yfe neighborhood of z = —7, the orbits lie on contours of constant H to 
0{e) for times 0(e~ 2). In the original three-dimensional map, the pendulum's librating orbits 
correspond to tubes aligned with the y-axis; two such orbits are shown in Fig. [3] (the red and 
purple orbits). The rotating orbits of (12) correspond to rotational invariant tori of (fsl) . The full 



width of this resonance in z is the width of the separatrix contour of H , appropriately scaled in e, 



(13) 



W(M,a) = ^\Ma\ . 
When ea = 0.005 as in Fie. [3] this gives w a = W(l, a) ~ 0.11. 




Figure 3: Phase space of M for /3 = 2, 7 = \{s/§ - 1) w 0.618, 5 = 2<5h « 0.0284, a = 6 = c = 1, and e = 0.005. 
Orbits trapped in four resonances are shown. The three blue-toned orbits correspond to the two (0, 1,0) resonances 
at z± = ±0.119. The green orbit is trapped in the (1, —1, 0) resonance at z c = —0.371, and the two red-toned orbits 
librate in the (1, 0, 0) resonance at z a = —7 = —0.618. Also shown are three rotational tori and a chaotic orbit (grey). 



A similar analysis can be done for the (0,1,0) resonance. If 5(3 > 0, the frequency map ([7]) 
crosses this resonance at two positions, z^ 1 , given in Tbl. [lj This resonance has rank one when 
Do = Q x (zb) = 7 + Zb Q, in which case the same expansion procedure can be applied for the 
variables (£,77, C) = (x,y,£^ 1 ^ 2 (z — z^)), since the fast angle is now x. If the fast frequency is 
Diophantine, Thm. pj applies and the averaged Hamiltonian is again (12) with new arguments: 
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H (y, £; Mf,, b). The effective mass Mb, given in Tbl. [TJ is real precisely when the resonance crossing 
is transverse: 5(3 > 0. In this case the resonance width is Wb = W(Mb, b) with the same function 



(13). Since the effective masses of the z^ resonances have opposite signs, the resonance with 
Mbb > is an island centered at y = 0, while that with Mbb < is phase shifted, and is centered 
at y = ^. The three blue-toned orbits shown in Fig. [3] are trapped in these resonances. For the 
parameters of this figure the resonances are centered at = ±0.119 and have widths Wb ~ 0.16. 
Since the separation between the resonances exceeds the sum of their half-widths, the Chirikov 
overlap criterion leads to the expectation that there are rotational invariant tori between the two 
resonances [LL92] . Indeed, this is what we observe numerically; one such separating torus (in grey) 
is shown in Fig. [3} 

The frequency curve @ also crosses the (1,-1,0) resonance transversely at the two points 
z = zf given in Tbl. [l] provided that P(6 + i)>-\. This resonance can be treated as before upon 
defining new variables 

(^r ] ,C) = (y,x-y,e- 1/2 (z-z c )) , (14) 
aligned with the resonance. The map Q now becomes 

£' = £ + u;o + 0(^) , 

with ljq = 7 + z c and M c given in Tbl.[l] Whenwo € V(c,fi), averaging over £, according to Thm.[5j 
now eliminates the terms proportional to a and b in g. The resulting system is again approximated 
by the time y/e map of ( 12 ): H(rj, (; M c , c). In Fig. [3] only the resonance at z~ « —0.371 is visible; 



it has a width w r . « 0.072. 



c 



w b 


+ w b 


M 


~ z b\ 


w t 


+ w~ 




- Zc\ 


w c 


+ w a 


2\zc 


- Z a \ 



In the phase portrait of Fig. [3] the resonances are all separated by rotational invariant tori. 
Indeed the Chirikov overlap parameters between neighboring, forced resonances, 

wt + w7 
s b +b- = b zr: « 0.66 , 

2 K - z b I 

s b -c- = ^ ~ 0.47 , 

0.37 , 

are less than the phenomenological threshold value s = | for "global chaos" |LL92| IMei07j , so we 
expect — and observe — that rotational invariant tori separate these resonances. Of course, there are 
also many chaotic trajectories for the parameters of Fig. [3j For example, the (1,0,0) resonance 
is surrounded by a large chaotic layer (grey in the figure) because of overlap with a (1,-1,-1) 
resonance at z* = —0.691 that is not shown in the figure. This resonance and the (1, 0, 0) resonance 
have an overlap parameter s = 1.2, and are not separated by any rotational tori. 

As 5 is varied the resonance locations and widths change as shown in Fig. |4j The averaged 
Hamiltonian description provides a qualitative description near each transversal crossing provided 
that the resonance widths are small enough to avoid overlap. Of course, even in this case the region 
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near the separatrix of the averaged Hamiltonian is replaced by a chaotic zone for the full map. As 
we will see in || a similar structure governs the dynamics near any rank-one resonance whenever 
Q(z) crosses the resonance transversely. 




Figure 4: Resonance centers z* (curves) and widths (grey regions) as a function of S, for the parameters of Fig. [3] 

However, when the resonance intersection is nontransversal, the dynamics is radically different. 
The curves z*{5) for the (1,0,0) and (1,-1,0) resonances in Fig. [4] are parabolas with turning 
points, z*{5t) = zt that signal nontransversal crossings. Tangency with an (m, n) resonance 
occurs when 

m ■ £1(zt) = n and m ■ DQ(zt) = . (15) 

Since a frequency curve £l(z) will generically intersect a resonance line for some z*, tangency is 
a codimension-one phenomena — it will generically occur if the frequency map £1 depends upon 
a single parameter. For the example Q we can think of 5 as the parameter that unfolds these 
tangencies. For the model ([8]), the (1,0, n) resonance is never tangent but tangencies do occur for 
both the (0, 1, n) and (1, —l,n) resonances where their effective masses diverge. The values of 5 for 
these tangencies are given in Tbl. [TJ 

The frequency curve is tangent to the (0, 1,0) resonance when 5 = 5t = at z = zt = 0. If we 
assume that Q, x {zt) = 7 G P(c, /i), then the angle x rapidly rotates, and for z « zt and 5 w 5t, 
the map Q is approximated by the averaged map 

y' = y — 5 + 0z' 2 , , , 

, (16) 
z = z — ebsm(2iry) . 

Again the averaged system is a two-dimensional, symplectic map. This system should approximate 
the original dynamics when z and 5 are near zt = St = 0. When j3 = 0(1), a maximal balance of 
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the terms in this map occurs when 



5 = 5t + e 3 A . 



2 

In this case, (16) is approximately the time £3 map of the Hamiltonian 



For |A| > \Ar\ where 



H R (y, C; 0, A, b) = ^ - A( - A cos(2vry) . (18) 

When A/3 > this Hamiltonian has four equilibria, at y = 0, | and £ = ^\J~^- When fib > the 
two equilibria (0, — \J~^) and (3; are saddles, and the remaining two are centers, see Fig. J^J 

9 

16tt 2 

there are two independent island chains. The separatrices of these islands correspond to the energies 

3 3 

Ee T = db — ?=(Ap — A2). These coincide at A = As, where there is a "reconnection bifurcation", 

see Fig. pj Note that reconnection occurs when sgn(A^j) = sgn(/3) at an unsealed value 5r = 
5t + s^Ar. The elliptic and saddle equilibria on the lines y = and y = \ move together as A 
decreases until they are destroyed in saddle-center bifurcations at A = 0. 

This structure can be most easily seen in the three-dimensional map by using slices. A slice is 
an approximate, two-dimensional phase portrait obtained by plotting orbits only when they land 
within a thin slab. In Fig. [6j we show two such slices, \x\ < 0.001 on the left and \y\ < 0.001 on the 
right — thus on average to obtain one point in either slice the map must be iterated ~ 500 times. 
These slices show the same orbits. The fast dynamics of the (0, l,n) and (1, — l,n) resonances are 
transverse to the left slice, and so the islands described by the slow dynamics can be seen. In 
particular the two (0, 1, 0) resonances at z^ are well described by the upper portrait of Fig. [5] since 
St = 0, and 5 = 25r at these parameter values. There are also three (1, —1, n) resonances indicated 



in the figure with their widths computed from (13). The (1,0, n) resonances do not appear in this 



slice since their librating orbits move approximately parallel to the slice; however, the chaotic orbit 
above the (1, —1, —1) resonance is due to its overlap with (1,0,0). Two (1,0, n) island do appear 
in the right slice; the diagonal (1, — l,n) resonances also appear again. 

Similar pairs of slices of the three-dimensional map are shown in Fig.JTJand Fig.|8]for 5 = 5r and 
\$r, respectively. Again, the averaged Hamiltonian portraits from Fig. 5 are closely mimicked near 
the (0,1,0) resonance in these figures except that the separatrices become chaotic. In particular, 
note that the invariant tori trapped between the resonances in Fig. [8] are "meandering" — they 
are not graphs over the unperturbed angles. 

Tangency similarly occurs for the (1,-1,0) resonance when 5 = 5t as given in Tbl. [l] when 
z = zt = jg, where the effective mass M c is infinite. As before, we must use the aligned angles 



(£, rj) of (14) to construct the averaged map, now scaling z and 8 as in (17). The result is exactly 
( [181 ) w ith parameters Hr(t], (, — /3, — A, c). Thus the same reconnection scenario applies to this 
resonance near 4/3(5 + 7) = —1. A three-dimensional phase portrait near the tangency is shown 
in Fig. |9j For these parameters, 5t ~ —0.743, and <5r ~ —0.729. The figure corresponds to 
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Figure 5: Contours of the Hamiltonian ( 18 1 for three values of A when /3b > 



0.21A^j. The chaotic layer around the (1,-1,0) resonances overlaps with that around the 
1) resonance, as indicated in Fig. |3j 
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Figure 6: Two dimensional slices \x\ < 0.001 (left pane) and \y\ < 0.001 of the map |8| for the parameters of Fig. [5] 



The arrows indicate the widths computed using (131. The a;-slice of the (0, 1,0) resonances resembles the top pane 
of Fig. [| 




Figure 7: Two dimensional slices of the map |8| with the parameters of Fig. [6] except S = Sr w 0.0142. The a;-slice 
of the (0, 1, 0) resonances resembles the middle pane of Fig. [| 
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Figure 8: Two dimensional slices of the map qSp with the parameters of Fig. [6] except 5 = |Aj? « 0.0071. The 
x-slice of the (0, 1, 0) resonances resembles the bottom pane of Fig. pM 




Figure 9: Orbits of the map |s]) for j3 = 2, 7 = §(\/5 - 1), 5 = -0.74, e = 0.001, a = 4 and b = c = 1. Shown 
are orbits in the two, nearly tangent (1,-1,0) resonances (shades of red and light blue) and in the nearby (1,0,1) 
resonance (green and dark blue). A chaotic orbit (grey) wanders among all three resonances and an invariant torus 
(yellow) bounds the trajectories below. 
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4 The General Rank- One Resonance 



In this section we consider more generally the case of a rank-one resonance of a volume-preserving 
map f e : T d x R — > T d x R that is an e <C 1 perturbation of the integrable map pb: 



x' = x + Q(z') + eX(x, z; e) , 
z = z + eZ(x, z; e) . 



(19) 



We suppose that (19) is volume preserving, so that its Jacobian has determinant \Df e \ = 1. In 
order to apply the averaging theory (see Thm. [5]), we will assume that X, Z, and Q are analytic 
functions of all their arguments. We will also need a Diophantine condition on the fast frequency 
vector at resonance, see below. 

When a point z* corresponds to a rank-one resonance, the frequency vector Q* = £l(z*) obeys 
a resonance condition 

m ■ 0* = n 



for exactly one coprime vector (m,n) £ Z d x Z \ {0}. Our goal is to show that (19) reduces 
essentially to ^ (with a general function g) near z = z* . 

As we discussed in the frequency curve will generically intersect a resonant plane transversely. 
However, under variation of a parameter, the intersection may become tangent; we assume that 
any tangency is quadratic. These two cases can be treated simultaneously by expanding about the 
resonant torus, setting z = z* + 0(e p ), and choosing the exponent p to be \ or \i respectively. 
In order to unfold the codimension-one tangency, it is convenient to expand about a torus slightly 
shifted away from the resonance. Thus instead of expanding about z* , we expand about a shifted 
point z, 

z = z + e p C, U<p<\. (20) 
We also assume that z = z* + 0{e 2p ) to define the frequency shift 5 by 

m ■ Q(z) = n- e 2p 5 . (21) 

Substituting these expansions into (Il9|) gives 



x 



' = x + n(z) + e p Dn(z*)C + U 2p D 2 n(z*)(' 2 + 0(e, e 3p ) , 



(22) 



C' = ( + e 1 - p Z(x,z*;0) + O(e) . 

Though the components of the vector (m, n) are coprime by assumption, the components of 
m need not be coprime, see App. |Aj Let k = gcd(m) be their greatest common divisor so that 
m = m/k is a coprime, integer vector. As discussed in App. [Aj there exists a basis for the lattice 
Z d with m as one basis vector; equivalently, there exists a matrix M G SL(d, Z) whose last row is 
rh T . Thus we can define new coordinates (£, rj) £ T d_1 x T as 



= Mx 

to align the last angle r/ with the resonance. 



Mx 
rh T x , 
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Since Z(x,z;e) is analytic, it has a Fourier series with analytic coefficients Zj(z;e). Conse- 
quently the perturbation Z(x, z*; 0) can be written in terms of the new angles as 



g(£, V )^Z(M- 1 (e,ri) T ,z*;0)= £ Zj< 



2ni(Z T ,r,)M- T j 



E 



9i e 



2rf-(£ T \r)) 



so that g has the Fourier coefficients gi = Z m ti, I € Z rf . 
Since z = z* + 0(e 2p ) by (21), we obtain 



kMn(z) 




where O* = kMVt* . In the new coordinates, the map (22) becomes 

e' = e+i(^+^Ac')+o( £)£ 2p ), 

r/ = r, + i (n + e p a(' + e 2p (P(' 2 - 6)) + 0(s, e 3p ) 



where 



1 



A = kMDU(z*) , a = m-Dn(z*), (5 = -m ■ D z Q(z 



(23) 



(24) 



Iterating this map k times to 0{e) will eliminate the rational § in the r\ equation, since this 
angle is taken "mod 1": 

Z k = Z + n* +£ P\( k + O(e,e 2 P) , 

Vk = Vo + e p aC k + e 2p (PC 2 k ~ 6) + e 3p ) , (25) 
a = Co + e 1 ~ p 5(6,w) + C»(£) , 

where 



fe-i 

E 



5 (£ + i^*,r/+ fn) 



This near-resonance map is partitioned into fast, £, and slow, (??, C)> variables. 

Note that (25) is essentially the same as ([8]), though with a more general force g. That is, 
if one neglects the higher-order terms in (25) and then formally undoes the scaling, by setting 
(x, y, z) = (£, i], e p C), the result is a d-dimensional extension of ^ (with an additional parameter A 
that can be also removed by scaling when it is nonzero). This reinforces calling Q the "standard 
volume-preserving" map. 

Since Q* satisfies precisely one resonance condition and the columns of M T are independent of 
rh, whenever I G Z is nonzero I • Q,* = kl ■ Mtt* = k{M T l) ■ n* $ Z. In particular none of the 
components of this transformed frequency can be rational: the angles £ are "rapidly varying." If, 
in addition, Cl* is Diophantine, Thm. [H] implies that, on the time scale t = 0{e p ~ l ) S> 1, the fast £ 
angles can be averaged away, reducing the system (25) to the two-dimensional map 



rj k = m + e p a( k + e 2p (f3( 2 - S) + 0(e, e 3p ) 
(k = (o + e 1 - p g( V o) + 0(e), 



(26) 
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where 

fe— l 

contains only the Fourier coefficients of g that are multiples of (0, k). Thus 

9~{v + l) = g{v) ■ 

Note that (26) is area preserving. In jj5]we will show that this property holds more generally, and 
that the degree of approximation between the averaged and the original system can be increased 
with additional coordinate transformations. 

To lowest order, the averaged map can be related to a Hamiltonian flow. For this to be possible 



we assume that go = 0, i.e., that g has zero average. Even though the derivation of (26) is valid 
when this is not true, the resulting area-preserving map would then have nonvanishing net flux and 
consequently, no invariant tori; it also would not be related to the flow of a Hamiltonian. For the 
rest of this section we assume that the average of Z is zero, as we also did in <|3j 

There are two distinguished limits corresponding to transversality or tangency of the frequency 



curve with the resonant plane; these can be treated by selecting the exponent p. If a, (24), is 
nonzero, the frequency curve crosses the resonance transversely (since the vector m is perpen- 
dicular to the resonant plane, and DO, is tangent to the frequency curve). In this case we set p = \ 



and the averaged map ( 26 ) is approximately the time \fe flow map of the Hamiltonian 



with DV(rj) = —g(rj). This is similar to the standard pendulum (12) and gives a resonance of size 
yje about z* in the original coordinates, as discussed in §[3j 



On the other hand, if a = 0, the frequency curve is tangent to the resonance plane. If (3, (24) 

is nonzero, then its curvature at the tangency is nonzero. In this case we set P = \ and (26) limits 

2 " — 

on the time es flow map of the Hamiltonian 

H R (r ] ,C) = ^C 3 -5C + V(r ] ) . 



This is similar to the standard reconnecting Hamiltonian (18). 

We have already studied in detail the two cases a/0 and a = 0, (3 ^ and their corresponding 
averaged Hamiltonians in Consequently, these examples describe the typical structure of the 
dynamics near a rank-one resonance for the cases of transverse and quadratically tangent reso- 
nances. Of course, in ^3] we only considered the dynamics with a specific, simple choice of g. What 
is surprising is that the numerical comparisons in ^3] indicate that the averaging theory (which is 
only valid for time of order e p_1 ) nevertheless gives a qualitative understanding (to the level of 
resonance structure, widths and the onset of chaos through overlap) of the long time dynamics. 
There will, of course, be many features that this simple theory misses. 
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5 Averaging for Volume-Preserving Maps with Nearly Constant 
Frequency 

When a dynamical system can be separated into slow and fast components, averaging can often be 
used to describe the slow dynamics on a long time scale [LM881 [SVMn7] . In this section we prove 



an averaging theorem relevant to the system (25) for which the fast dynamics, the £ variables 



corresponds to translations on a torus with nearly constant frequency. In fact, the theory applies 
to more general maps with any number of fast and slow variables. To reinforce this, we denote 
the fast angles by 9, and the slow variables ((ry, Q in the previous section) by J — to indicate their 



similarity to "actions". The slow- fast system (in particular also the map (25)) can then be written 

a' = e + i0o + KF w (9,J;K) , 

(27) 



J' = J + kG W {6, J;k) , 



on T x M where M CW for some j. Here k represents the small parameter e p of (|25|) . 

We will show, in addition, that if the original system is volume preserving, then the coordinate 
transformations that achieve the (higher-order) averaging can be chosen to be volume preserving as 
well. The implication is that the truncated, averaged system will have a symmetry (i.e., translation 
along the fast angles) and thus the decoupled, slow system will again be volume preserving. In 



our particular example this shows that the map (26), in appropriately transformed slow variables 
(77, (), is area preserving. 

The goal of the averaging is to obtain a description of the slow dynamics on the time scale 
0(fv -1 ). One hopes that if k <C 1 and ujq is sufficiently incommensurate, then the slow dynamics is 
approximated by 

J' = J + kG {1 \J) , (28) 

where 

G(J) = [ G(8,J;0)d9 . (29) 
This first order averaged equation should give dynamics close to that of the full system in the sense 



that for any orbit (6t, Jt) of (27), then provided that Jo = Jo we should have 

\J t -J t \=0(K) 

for < t < ^ for any fixed T > 0. 

In most cases averaging theory is formulated for ODEs, though there are some results specifically 
for maps [Kif03, DEV04, Bra09]. The analyses of Kifer and Brannstrom apply to systems in which 



the slow and fast dynamics are more strongly coupled than (27) and for which the slow dynamics 
need not correspond to toral translations. In [Bra09] the angle dependence of G is assumed to occur 
only at 0(k 2 ). These theorems apply only when the constant vector loq is replaced by a function 
f2(J) that satisfies a nondegeneracy, or twist condition (Kolmogorov nondegeneracy) . In either 
case, they must deal with the difficulty that the ergodicity of the fast angle dynamics can depend, 
at lowest order, on the slow variables; consequently in order to apply these theorems, a nontrivial, 
"good-ergodization" property must be verified. The results are weak in the sense that the averaged 
and true, slow dynamics remain close only in the sense of measure. The reason we can achieve 
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better results is that we are treating the simple case where the frequency of the unperturbed system 
is constant and Diophantine. Because of this we are able to achieve arbitrary high order accuracy 
for every initial condition. 

A standard implementation of averaging begins with the construction of a coordinate change to 
push the coupling to the fast angles off to higher order |SVM07| . Being ultimately only interested 
in the slow dynamics, we could choose not to transform the angle dynamics. However, we will 
see that the truncated, slow map is volume-preserving if the fast angle variables are transformed 
as well. When loq is Diophantine, Q, and the map is analytic, this transformation can be done 
recursively to any order. We start by assuming that the normalizing transformation has been done 
up to 0(n n ~ 1 ), but the result also applies to (27) directly by setting n = 1. 



Lemma 1 (Decoupling Transformation). Let f be the real analytic map on T l x M 

n-1 



n-1 



J' = J+Y1 Ki G U) (J) + K n G^(6, J; K ) , 

i=i 



(30) 
(31) 



such that ojq G T>(c, ji) for some c, fx > 0. Then there is a kq > such that for all \k\ < kq there 
exists a diffeomorphism h^ n \9,J) = (inversely) defined by 



(8, J) = (V + « n rM(^ I), I + n n S^ (V, /)) , 
that conjugates the dynamics to the analytic map 

n 

= $ + oj + KiFU) ( J ) + (^, I; , 

n 

I' = I + J2 K j G U) (I) + K n+1 G {n+1) (if>, I; k) , 

3=1 



(32) 



(33) 



where G n and i* 1 ( n ) are the 



averages 



f Q( n ) and F^ n \ respectively, as defined by (29). 



Proof. The elimination of ^-dependent terms to 0(k h ) proceeds exactly the same for the 9 and J 
equations. We present the argument for (31) — replacing (G,S) by (F,T) would give the argument 
for (30) (the affine term, uq, simply carries through). 

To eliminate ^-dependent terms substitute (32) into (31) and expand to n th order to obtain the 
homological equation 



S^(^ + u ,I)-S^(^,I) = G ( - n \^,I;0)-G^(I) . 
Formally, this can be solved by Fourier transformation, setting 



(34) 



I) = S£° {I)e 2 ^ , G< n > (i/>, I;0)=J2 (I) e 2 ^ , 



kei 1 
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and noting that G^\l) = G^ n \l) to obtain 



(35) 



The function S^ n) (I) = S^ n \l) 



y± ) — u • ■ \± j can be freely chosen, for example, to be 0. 



Since loq is incommensurate, the denominator of (35) never vanishes: each Sk(I) is an analytic 



function of /. Moreover, when ujq is Diophantine and G is analytic in 9 then the formal Fourier 
series for S converges to an analytic function, see for example [Mos66]. 
The new map for / is 



I' = J' -K n S (n \lp',l') 



n-1 



I + Y, ^ j G u) (J) + K n [G^ n \e, J; k) + S( n \il>, I) - S^(iP', I')} 



(36) 



Comparing with (33), gives an implicit definition for G^ n+1 ^: 



G^ n+1 \i>,I;K) = - 



,j-n 



n-1 



J) - G (j) (/)} + G (n) (9, J; k) - G (n) (I) 



+ 



(37) 



where J and 6 are just short-hand notation for the right hand side of (|32j). A similar equation 
is obtained for F( n+1 ) and these equations are implicit and coupled because is a function 

of (G (n+1) ,F (n+1) ). The homological equation Q implies that G (n+1) is well-defined as k -> 0, 
because it forces the O(kP) term inside the square bracket in (37) to vanish. The leading order 
term at k = thus is 

G {n+1) {6, 1- 0) = DjG , W(J)S'W(0, /) + D K G^ n \e, /; 0) 

- D e S^ {e + Lo , 1)PW (/) _ Dl S^ (9 + ojo, T)QM (I) . 



Moreover, since the Jacobian of the right hand side of (37) and the analogous equation for 

with respect to (G^ n+1 \ F^ n+1 ^) vanishes at k = 0, the implicit function theorem implies that there 

is a kq > such that for all |k| < kq, the functions G^ n+1 \ F( n+1 ) exist and are analytic. 



Finally, since the transformation (32) is a near-identity transformation, kq can be selected so 



that it is a diffeomorphism. Thus the dynamics of (ift, I) and (9, J) are diffeomorphically conjugate. 

□ 



We now show that (32) can be replaced to by a volume-preserving transformation when / 



is itself volume preserving. In order to do this we interpret (T< n ), S< n )) as a vector field V 
(T^d^p,S^di). This vector field generates a transformation e~ v (if),I) that is volume preserving 
whenever V is divergence free, V • V = 0. To show that this is possible, we first prove two 



simple lemmas. The first exploits the assumption that the Jacobian of (27) has the form Df = 
I + kB\ + 0(k 2 ). If / is volume preserving, then 1 = \Df \ = 1 + ntr B\ + 0(k 2 ), which implies 
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tr B\ = 0. Note that higher order terms do not necessarily have zero trace, however. Put differently, 
when tr B = 0, then | exp B\ = 1, but only the linear term in exp B = 1 + B + \B 2 + . . . necessarily 
has vanishing trace. More generally: 

Lemma 2. Consider a convergent sequence of nonsingular, square matrices M n = Ylk=o K Bk such 
that the determinant \Moo\ = 1. Then for any finite n we have tr B n = C n (B n _\, . . . , B\) where 
C n is some smooth function of the matrices B^ with indices k < n. 

Proof. Firstly, writing M m = M n _i + 0(n n ) and taking the determinant of both sides gives 1 = 
|M n _i||7 + 0{K n )\ = |M„_i| + 0(k h ), which implies 

|M n _!| = 1 + 0(k u ) = 1 - K n C n (B n ^,B n . 2 , . . . , Bi) + 0{K n+1 ) . 

Secondly, observing that M n = M n _i + K n B n and taking determinants of this identity gives \M n \ = 
\M n _i\\I + K n M~\B n \. Using the previous estimates for the determinants of M n and M n _i then 
implies 

1 + 0(K n+1 ) = (1 - K n C n + 0(K n+1 ))(l + n n tr B n + 0(K n+1 )) , 

where in the last term we have used that M n _i = I + O(k) so that M~ j = I + O(k) and hence 
tr(M~\B n ) = ti B n + O(k). Finally, expanding the right hand side and equating terms of order 
0(n n ) in this equation gives tr B n = C n . □ 

Next we establish a connection between the divergence of a periodic vector field and the diver- 
gence of its average. 

Lemma 3. Let V be a vector field onT 1 x M in coordinates (9, J) and denote by a bar the average 
over all angles 6. IfX7-V = C(J) where C is a function of J only, then V • V = V • V. 

Proof. First we show that in general V • V = V • V. Write V = (F, G) where F : T l x M -> T l and 
G : T l x M — > M so that V ■ V = ti DgF + tvDjG. Since F is periodic in 9, the average of any 
0-derivative vanishes. Hence V • V = tr DjG = tr DjG = V • V since the operations of averaging 
over 9 and differentiation with respect to J commute. By assumption we have V • V = C{J) 
and averaging this identity gives V • V = C(J) = C(J). Finally, replacing V • V by V • V and 
eliminating C(J) gives the result. □ 

With these two lemmas we can now prove: 

Lemma 4 (Volume-Preserving Decoupling Transformation). If the map f is volume-preserving 
then the decoupling transformation (32) can be extended to a volume-preserving map. 



Proof. The volume-preserving extension of the near-identity transformation (32) is constructed by 
interpreting V = {T^d f , S^d T ) as a vector field and using its exponential as the transformation 
operator. What we need to show is that if the map / is volume-preserving then V can be chosen 
to be divergence free, so that the resulting exponential is volume-preserving. 

The first observation is that the homological operator (34) considered as acting on vector fields, 
preserves the subspace of divergence free vector fields and its complement. Introducing the unper- 
turbed map h 0J0 (9, J) = (9 + lhq, J) the homological operator can be written as 

L^V = {Dh wo y l V oh-V= (h* - id)V . 
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Now if V • V = 0, then this also true for the image L aJo (V) because h WQ is volume-preserving. 
Moreover, since L Wo is semi-simple, the complement of any invariant subspace is preserved as well. 

To complete the proof, we need to show that the right hand side of the homological equation, 

is a divergence free vector field, so that the generator of 
the transformation (T^ n \ S^) can be chosen to be divergence free as well. 

To apply Lem. |jwe identify the Jacobian of ([30|) and ^ with M n , i.e., Bj = D{F^\G^) 
for j = l,...,n- 1 and set B n = D(F^ n \i/j, I; 0), G^^, I; 0). The result is that tr B n = V • 
(F( n )(^, I; 0), G( n )(^, I, 0)) = C n (5 n _i, . . . , Si) = C n (I) since all lower order terms have already 
been averaged. This sets the stage for Lem. [3J which implies that V • (F^ n \tp, I; 0), G^ n \tp, I;0)) = 
V-(fW,GW), so that the divergence of the difference (fW-FW, GW-GW) vanishes. Therefore, 
the right hand side of the homological equation is divergence free, and hence the generator (T, S) 
can also be chosen to be divergence free at each stage. 



Finally, if we replace the transformation (32) by (e 



I), we obtain the result. □ 



Now we return to the general (not necessarily volume-preserving) case. We will show that, 
roughly speaking, the projection ir(t/j, I) = I of the slow dynamics of (27) onto M is conjugate, to 
0(n n ), to the averaged map 



(38) 



for times of order k^ 1 for any n > 1. More precisely: 



Theorem 5 (re Order Averaging). Suppose that f is an analytic map of the form (21) and uq is 
Diophantine. Let h = o o . . . h^ 1 ) be the composition of the sequence of transformations, 

(iP,I) = h^(0,J) , 

given by the inverse of (32), which exists for \k\ < kq for some kq > by Lem. [7| If given a T > 0, 
the orbit It of (38) is in some compact subset K C M for < t < — , then 



\wh(6 t ,J t )-I t \=0(K n ) , 
where (6t, Jt) is any or&i£ o// with initial conditions satisfing Iq = irh(8o, Jq) 



(39) 



Proof. By Lem.JT] h conjugates (27) to (33). Thus 7rh(9t, Jt) = It, an orbit of the second component 
of g. Since is analytic on T_ z x K and if is compact, it is bounded, say |G( n+1 )| < C. 

Similarly, on K the function X^j=i ^G^(I) is Lipschitz, say with constant kL. Consequently, 



lit 



.n+l 



U\ <K\I t - 1 -It^ l \+2CK n 

and the discrete Gronwall lemma, e.g. [[Arg00|, implies that 

\It ~ h\ < 2CK n+1 t(l + kL) 1 - 1 < K n (2CTe LT ) 

for all < t < -. 



□ 
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Remark 1. It may seem difficult to find an appropriate set K that contains the orbits of interest 
over time T/k in the limit of vanishing k. For this it is useful to notice that for sufficiently small k 
the orbits of It are well approximated by the flow of the divergence free vector field for time 
T. If this time-T flow is bounded, then the region of M covered by the flow is a good candidate 
for the set K. 

In Thm.[5j the map was normalized up to 0(k u ) with remaining correction terms (_F(" +1 ) , G^ n+1 ^) 
of order 0(n n+1 ). In practice it is easy to compute the averages (_F(n+l) ) G( n + 1 )), which gives the 



next order term of the averaged map ( 38 ) whose orbits we denote by II . It turns out that it is 



not necessary to compute the next order transformation /j( n+1 ) to get an error of order K n+1 . 



Corollary 6. Denote by I^ n+l ^ be an orbit of (38) with n replaced by n + 1. With the other 



assumptions as in Thm. [5] we then have 

\TTh(e t ,Jt)-ii n+1) \ = o(K n+1 ) . 

Proof. If we apply Thm. [5] at order n + 1 the statement is 

|7r(/i( n+1 ) o h)(9 t , Jt) ~ I (n+1) \ = 0{n n+1 ) . 
Observing that h^ n+1 ^ = id + 0(k u+1 ) and expanding the right hand side gives the result. □ 
Combining nth order averaging, Thm. [5j and volume-preserving decoupling, Lem. [4j we find: 



Corollary 7 (Volume-Preserving Slow Map). If the map (21) is volume-preserving, then the trans- 
formation of Lem. [7] can be done so that the truncated slow map (38) is volume-preserving on M 
up to order n. 



Proof. After decoupling and truncating at order n, the transformed map (33) is independent of ij; 
(except for its identity part). Consequently, the Jacobian of the truncated map is upper block- 
diagonal, and the upper left block is the identity. This results from the translation symmetry of 



(38) along ijj, and its semidirect product form. Thus the determinant is equal to the determinant 



of the lower right block, which is the Jacobian of the slow map (38) alone. As in Lem. |4j the 
decoupling transformation h can be chosen to be volume preserving, and so the determinant of the 
Jacobian of the untruncated map is equal to 1. When terms of order higher than n are neglected 
the determinant is 1 + 0(K n+1 ). □ 

Remark 2. One may choose to leave the fast dynamics untransformed, setting 9 = i[> in each 
decoupling step. This would still decouple the slow dynamics from the fast dynamics, and Thm. [5] 
still holds. However, the reduced map would not in general be volume preserving. 

We now give some details of this averaging procedure applied to the example ^ with perturba- 
tion Q near the (1, 0, n) resonance, where z* = n — 7. In this case the function g is a trigonometric 



polynomial. After blowing-up near the (1, 0, n) resonance, the map ( 10 ) has fast variable 9 = £ with 
wo = /3(n — 7) 2 — 5, slow variables J = (n, £), and small parameter k = ^fe. The fast dependence 
occurs in every component of the map; however, in the £ and the n equations, it only appears at 



quadratic order in k. As a consequence, the first order decoupling transformation (32) only changes 
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Figure 10: Projections often orbits of ^ with initial conditions x = J\ = —0.49 and £ = J2 = 0.02 
and covering y = in equidistant steps. Parameters are chosen as in Fig. [3] except 5 = 0.7. The 
panes in a row show irh(6t, Jt) for h = id, h = h = o faW from left to right, respectively. 
For the top row e = k 2 = 0.0002, t = 0, . . . , 177, while for the bottom row e = k 2 = 0.0001, 
t = 0, . . . , 250. The averaged orbit, of ( 11 ), is essentially identical to the lower right plot. 



Since g is already given as a Fourier series it is easy to compute the transformation: 

1 



0,0, 



2 sin(-7ra;o) 



(6cos(27r£ — ttujq) — ccos(27r(^ — 77) — ttujq) 



(40) 



Note the occurrence of the "small-denominator" sin(7ra;o) that arises from (35). Since the O(k) 
terms in are independent of £, S''- 1 ) is independent of £, and the resulting transformation is 
exactly volume-preserving. In the new variables the dependence on the fast angle £ appears at 
0{n 2 ) in every equation. Averaging these equations over £ gives the standard map (11). 

The next order transformation depends on all variables, and has all components non-zero: 




1 



4 sin ttujq 



2/3(n — 7)(5sin2-7r£ + csin2ir(r] — £)) 
b sin 27r£ + c sin 27r(ry — £) 
\2vrC(26/3(7 - n) cos 2vr£ - c(l + 2/3(7 ~ n)) cos 2^(1? - £)) / 



The corresponding transformation removes the £ dependence up to 0(k 2 ). Averaging the trans- 
formed equations at C(k 3 ) gives a correction to the standard map that has the sole effect of changing 
the coefficient of sin2-/rx from —na to 



bcir 



—na 



-na — k 



4 sin ttloq 

At the next order, additional terms of order k 4 in the truncated decoupled map are introduced and 
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the final averaged map becomes 



I[ = h + Kl' 2 + K 4 X (h) , 

I' 2 = I2 — Kasin27r/i — k 4 /2-D , 

where x(-^l) = 7i + a l sin27r/i + j3\ cos27rii. 

The orbits of the full map in the transformed coordinates are shown in Fig. 10 showing that 
as the order of the transformation is increased, or the size of the perturbation is decreased, the 
fluctuations due to the fast dynamics become smaller. Log-log plots of error vs. k clearly verify 
that the scaling of the deviation from the averaged orbits has exponents 2, 3, 4 as promised by 
Cor.! 



6 Double Resonance 

For the one-action map ([!]) in three dimensions, the image of the frequency map is a two-dimensional 
frequency vector oj. If its resonance has rank two, then oj must be rational, say, 

u= (- - 
V d\ ' d 2 



Note that orbits of the integrable torus map (43) with this rational frequency are periodic with 



period lcm(di,d 2 ). In this case, the near-integrable map (19) has two slow angles, and as we will 
see, the map cannot be reduced in dimension by averaging. 

To obtain the normal form, we first address the problem of finding an appropriate basis for the 
module C(oj). Note that, without loss of generality, we can assume that gcd(k,di) = 1; however, 
the denominators d\ and d 2 need not be coprime. When gcd(di,d 2 ) = 1 the only way to achieve 
m ■ oj G Z is to cancel the denominators di, implying that the resonance module has a basis 
mi = (di,0) and m 2 = (0,d 2 ). By contrast, when 

g = gcd(di,d 2 ) / 1 

the problem of finding a basis for C{oj) is more interesting. One basis vector is easy, either (dj.,0) 
or (0, d 2 ) suffices. However, together these two vectors do not form a basis, since there is additional 
cancellation. If for concreteness we choose mi = (di,0), then upon writing di = gdi, the equation 
to solve for the second vector m 2 = (p, q) becomes 

h h 

gm 2 ■ oj = —p + —q = gn, p,q,n G Z . 
d\ d,2 

This implies that p = d\i and q = d 2 j, so we arrive at 

hi + hj = gn, i,j G Z , 

where by construction gcd(li,g) = 1. To make m 2 independent of m\, we must choose j 7^ 0. We 
can always find a solution with j = 1, since the Diophantine equation 

gn — hi = l 2 
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is always solvable for integers i and n when gcd(g,l±) = 1 [HW79J. Thus a second basis vector 
is 77i2 = (idi,d2). Note that the determinant of the matrix of basis vectors is d\d,2 = gd\di = 
lcm(di, c^), which is the period of orbits with rotation vector co. 

Example. The module for 10 = (g, g) has g = 1, so we can take mi = (3,0) and 7712 = (0,2) and 
then 

£ = {(3i,2j),(i,j) eZ 2 } . 

Of course any unimodular transformation of mi and 771-2 also gives a basis; however, the selected 
basis is minimal in the sense that no other basis has smaller 1-norm. 



Example. The module for ui = (2; 4) has period 4 and g = 2. In this case the choice mi = (2, 0) 
and 771-2 = (0,4) misses some lattice points, for example (1,2). The process outlined above gives 
7712 = (1)2) and the same mi, which is the minimal basis in the 1-norm. 



Example. If ui = (3^,5) then the period is 36 and g = 3. If mi = (12,0), we must solve the 
equation i + j = 3n. For j = 1 a solution is i = 2, so we have a second vector 7712 = (8, 3). The 
minimal basis is (4,-3), (8,3). 

Once a basis for the module is known, the analysis of App. [A] can be used to choose a set 
of angles on the torus that have diagonal dynamics. In particular, we can diagonalize the angle 

PKQ where P and Q are unimodular matrices, and, 
diag(A;i, /C2). Letting 77 



dynamics using the transformation ( 44 ) , M 



since we have only two angles in this case, K 
unperturbed angle dynamics on the resonant torus, $l(z*) = 

7/ = 77 + Q£T = 77 + Cj 



Qx, then reduces the 



n* of (19) to 



A 



k\ divides hi- 



where Cj\ = Note that this map has period /c2, since, as discussed in App. 

To find the map near a double resonance, we can now follow the same steps we used in £]4j 
expand (19) about the integrable torus using (20) with exponent p = \, set (5 = 0, and iterate the 



resulting system /c2 times. This gives 



7]' = 77 + ^feaC! 



+ 0(e) , 



(41) 



where g{rf) is obtained by iterating A>2 times along the unperturbed torus map 

k 2 -l 

3=0 

with g(rj) = Z(Q~ 1 7], z*; 0). Here we assume that the twist vector, a = k2QDQ(z*), is nonzero. 
This map is a near-identity map that is, to lowest order, the time yfe map of the flow of the 
incompressible vector field 

77 = a( , 

C = g(v) ■ 



(42) 
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Note that this system still has a three-dimensional phase space, and that the "force" g is a periodic 
function of two angles; however, these move along a fixed line on the torus. Since a is generi- 



cally irrational this line generically has irrational slope, and the system (42) is a "quasi-periodic 
pendulum". Indeed, introducing the lifted coordinate s E M. by r\ = as, then this system has the 
Hamiltonian 

where V is a quasi-periodic function of s with derivative DV(s) = —g(ais,a2s). This model also 
occurs as a normal form for symplectic maps near a twistless singularity in the frequency map, see 
|D1M06| . 

7 Conclusion 

We have studied resonances in near-integrable one-action, volume-preserving maps. Near a rank- 
one resonance, one component of the frequency is rational (in an appropriate basis), and when the 
remaining frequency components satisfy a Diophantine condition, there exists a volume-preserving 
coordinate transformation that decouples the fast angles from the slow angle and slow action to 
arbitrary order in a small parameter that represents the deviation from exact resonance. Truncation 
of the resulting map gives rise to a skew-product form in which the slow dynamics is given by a 



decoupled, area-preserving map; at lowest order, the result is (26). This map describes the dynamics 
of the full system on a finite, but long time-scale, as shown by Thm. [5] and Cor. [7j 

There are two typical cases. If the image of the frequency map crosses the resonance transversely, 



then the slow map is essentially Chirikov's standard twist map ( 11 ). On the other hand, if fi(#) has 



a quadratic tangency with the resonance, then the slow map is essentially the standard, nontwist 



map ( 16 ) 



These two cases both arise naturally from the three-dimensional map (|8j), which we propose to 
be the volume-preserving analogue of Chirikov's standard, area-preserving map. Indeed, expanding 
a general one-action, volume-preserving map, near a rank-one resonance gives rise to this standard- 



form, (25). Moreover, the standard map (|8j) itself yields this form when it is expanded about any 
rank-one resonance — it is recapitulated by the blow-up transformation. 

The averaging results show that, near resonance, the dynamics of the three-dimensional standard 
map is surprisingly similar to that of the projected, two-dimensional area-preserving maps. In 
particular we showed that Chirikov's resonance overlap criterion can be used to estimate parameters 
at which invariant tori breakup near a transverse resonance crossing. 

Of course, for times longer than the O^n^ 1 ) of Thm. [5] (and in particular when resonance overlap 
is significant), our simple averaging theory no longer applies and the dynamics must be treated as 
a intrinsically three-dimensional. For example, an understanding of the break-up of rotational tori 
requires a fully three-dimensional treatment [Meillj . 

Recall that in the symplectic setting, vanishing twist is related to local non-invertibility of the 
frequency ratio map [DIM06J . Such singularities could also be considered for the volume-preserving 
case. These singularities might be global — when the frequency ratio map is not one-to-one, or 
local — when the map is not smooth at some point. An example is the development of a cusp 
singularity, modeled by the unfolded cusp normal form: 

U(z) = (cj + z 2 ,5z + z 3 ) . 
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We plan to study such cusp singularities in a future paper. 

The results outlined above hold near a rank-one resonance. By contrast, a very different reduced 
map (41) is found near a rank-two (or double) resonance. It is, however, again a version of the 
standard map, but in this case has a quasiperiodic force. 



A Resonant, Integrable Torus Maps 

Consider the map 

6' = 9 + cj mod 1 (43) 

on T d . Suppose that uj GT d has a rank-r resonance, that is the resonance module ^ has dimension 
r = dim£(u;), where < r < d. Let mi G ...i = l...r be a basis for C{lo), that is, a set of 
independent vectors so that 

£(w) = span z (mi, . . . , m r ) . 
The module C is a sublattice of 'L d . It is said to be primitive if 

C = span R (mi, . . . , m r ) Pi 7L d , 

that is, every integer point in the hyperplane spanned by the vectors rrii is a point in C Let M be 
the r x d dimensional matrix 

M T = (mi,m2, . . . , m r ) . 

By assumption rank(M) = r. 

Following Lochak and Meunier [LM88J and the standard normal form construction for matrices 
MI5!)T. Chap. XL 7], there exist unimodular matrices P and Q such that 

M = PKQ . (44) 

Here K is the r x d diagonal matrix 

K = diag(fci, . . . , k r ) , 

where ki £ N and k{ divides kj for i < j. The integers ki are the invariants of C. They determine 
when the lattice is primitive: 

Lemma 8. A module is primitive if and only if all its invariants are one, ki = 1. 
Defining the r-dimensional vector n by rii = uj ■ mj, i = 1, . . . , r then 

n = Mu) = PKQu . 

If we define Co = Qu, then we have 

KCj = P~ l n = h ; 

thus Qji = j*- for i = 1, . . . , r, while obi for i > r are incommensurate. 
We use the d x d matrix Q to construct a new set of angles on T d 
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so that the map (43) transforms to tp' = tfi + a), or in components 



Hi 



i>'i = i>i + 77 , i = l,--.,r 
ip'k = i>k + & k , k = r + l,...,d . 
The orbits of the first r phases are periodic with period 

k r = lcm(/ji, . . . , k r ) , 



(45) 



and since by assumption, there are no additional resonance relations, the orbits of the last d — r 
components are dense on T d ~ r . Thus we may conclude that 

Lemma 9. When uj is r -resonant, i.e., r = dim£(a;), the orbit of each initial condition of the 
integrable torus map (43) is dense on a family ofk r , (45), disjoint, d — r- dimensional tori that are 
mapped into each other by iteration. 

Example. For d = 3, the frequency uo a = (\/5, 2\/5+ \, \/2) is 1-resonant: in order that m-oj = n, 
then ?7i3 = 0, mi + 2rri2 = 0, and ttt-2 S 2Z. Thus the module for this resonance is 

£K) = {(-4j,2j,0):j€Z}cZ 3 , 

which has dimension one with a basis vector m = (—4, 2, 0). This resonance is not primitive since 
the components of the minimal basis vector m are not coprime. Indeed the line mt contains the 
integer vector (—2, 1,0) that is not in C{uj a ). 

For this case, the normalization (|44h results in 



M 



Thus ft = Pn = 1, and 



4 2 = (1) 2 




PKQ . 



Cj = Quj 



so that we have Klu = ft as required. Note that the transformation Q is orientation reversing, 
det Q = — 1; however we could fix this by exchanging the last two columns. 

Example. The frequency cj& = (\/5, \/5, 2) has rank 2 with module 

£K) = {ji(l,-l,0)+j 2 (0,0,l) :ji ez) • 

Here 

\0 1 
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is a minimal basis for C. This module is primitive. 



Example. The frequency oj c = (2v2+ \, \ ) is 2-resonant. A basis is 



M 




with n = (l,2) rj 



The normalization (44) is 



P 




K 






(l 


2 


l\ 


Q = 








1 




V> 


1 


0/ 



so that 



( 1 \ 



n 



n 



The orbits of this torus map lie on two disjoint circles. 

Example. If w G I 3 has rank three, then its components are rational. Indeed, if rrii ■ to = ni for 
three independent rrii, then oj = M~ 1 n, where the 3x3 matrix M is nonsingular by assumption. 
Note, however, that this does not imply that the module is primitive. Instead, the torus decomposes 
into period k r orbits. 
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